Optimal potentials for temperature ratchets 
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In a spatially periodic temperature profile, directed transport of an overdamped Brownian particle 
can be induced along a periodic potential. With a load force applied to the particle, this setup can 
perforin as a heat engine. For a given load, the optimal potential maximizes the current and thus 
the power output of the heat engine. We calculate the optimal potential for different temperature 
profiles and show that in the limit of a periodic piecewise constant temperature profile alternating 
between two temperatures, the optimal potential leads to a divergent current. This divergence, 
being an effect of both the overdamped limit and the infinite temperature gradient at the interface, 
would be cut off in any real experiment. 

PACS numbers: 05.40.-a, 82.70.Dd 



I. INTRODUCTION 

Noise induced transport occurs in a broad variety of 
systems both in physics and biology [ij. In a micro- 
scopic system embedded in a thermal environment, di- 
rected transport generally requires two ingredients: (i) 
External sources driving the system out of equilibrium 
and (ii) broken spatial symmetry [l|. The main idea is 
to impede the thermal motion in one direction in order 
to obtain a net current in the other direction [3]. This 
mechanism is called the "ratchet effect". Applications 
range from particle sorting [3, \^ to modeling molecu- 
lar motors 0, [a, [y, 0| • Most studies on ratchet motors 
focus on a given potential landscape and a given driving 
scheme. For practical purposes, however, optimization of 
the driving mechanism with respect to a maximal current 
is an important issue. 

For discrete analogues of ratchet motors, where the 
potential landscape is characterized by only a few pa- 
rameters, such an optimization has been performed in 
models for microscopic heat engines [8[ and molecular 
motors 0. Paradoxical Parrondo games [10| which can 
be interpreted as discrete analogues of Brownian ratchets 
[ll| . have also been optimized recently |13 |. 

For continuous motors, where optimization requires 
variational calculus, there exist only few results. The 
optimization of driving schemes has been studied for 
time-dependently driven ratchet motors [ij, [lj| and for 
a Brownian heat engine cycling between two heat reser- 
voirs [131 • Potential landscapes have been optimized for 
the transport across membrane channels [iq . Maximiz- 
ing the current of flashing ratchets by using a feedback 
control strategy has been proposed recently [ij, [131 ■ 

In this paper, we focus on a continuous thermal 
ratchet where transport along a spatially varying time- 
independent potential is driven by a periodic spatial tem- 
perature profile [l^, [20, [21I, [24I . The recent generation of 
temperature gradients on small length scales [23, [2J, [23 
may render such molecular heat engines experimentally 
realizable. Cargoes driven by thermal gradients on a sub- 
nanometer scale have already been observed [26j . 

Thermodynamic efficiencies of such ratchet heat en- 



gines have been calculated in Rcfs. [27|, |28|, |29| . It has 
been argued that heat engines should generally be char- 
acterized by their performance at maximum power out- 
put [so, [3l| ■ Recent studies on ratchet heat engines have 
varied the load force for a constant potential in order to 
maximize the power output [32l . |33| . We take a comple- 
mentary approach and optimize the potential for a given 
load. The maximization of power output and particle 
current then is equivalent. 

So far, this task has been tackled numerically only for 
a one-parameter class of potentials for a given temper- 
ature [SJ]. To the best of our knowledge, there exists 
no systematic study on the optimal potential for such a 
thermal ratchet. Using variational calculus we determine 
the optimal potential which maximizes the current for a 
given temperature analytically up to a numerical root 
search. For a dichotomous temperature profile with in- 
finitely steep gradients as used in most previous studies 
[29l. |32|. |33|| . the maximal current diverges. This unphys- 
ical behaviour is an effect of the idealized assumptions 
of both an infinite temperature gradient at the interface 
and the overdamped dynamics. In any realistic system, 
temperature gradients will be finite which is sufficient to 
cut off the formal divergence. 



II. THE TEMPERATURE RATCHET AND ITS 
CURRENT 

We consider a Brownian particle of mass m moving 
in a periodic potential V{x) with V{x + L) = V{x) in a 
viscous fluid with friction coefficient 7. A constant load / 
is attached to the particle. The surrounding temperature 
is modeled by T{x) which has the same periodicity L 
as the potential. A special case is a piecewise constant 
temperature with a hot and a cold area, see Fig. [1] For 
a properly chosen potential the particle moves against 
the load on average. The thermal fiuctuations in the 
hot area can push the particle against the load over the 
barrier of the potential. As soon as the particle is in 
the cold area the probability of getting pushed back is 
smaller because of the weaker thermal fluctuations. In 



this way the particle drags the load and produces work 
effectively acting as a heat engine that works between two 
heat baths. Such a mechanism is not limited to piecewisc 
constant temperatures. 




FIG. 1: A temperature ratchet: a particle with a load / mov- 
ing in a potential V{x) in two piecewise constant temperature 
regions Th and Tc. 

The time evolution of the position of the particle x{t) 
is governed by the Langevin equation 



mi = -7i - V'{x) +f + g{x)£,{t), 



(1) 



where time derivatives are denoted by dots and space 
derivatives by primes. The stochastic term g{x)S,{t) mod- 
els the thermal noise from the environment. Its strength 



g{x) = y/kBT{x)^ 



(2) 



depends on the temperature profile T{x) resulting in mul- 
tiplicative noise. 

In the following we want to make two assumptions: 
the friction dominates the inertia and the noise is un- 
correlated Gaussian. When dealing with multiplicative 
noise, these two limiting procedures do not commute 
[35l . [Sq . The order of the limiting procedures determines 
the stochastic interpretation of the term g{x)^{t). If we 
first assume Gaussian white noise 

{at)m)^mt-t'), (3) 

and then go to the overdampcd equation 

-fx^~V'{x) + f + g{x)at) (4) 

we end up with a corresponding Fokker-Planck equation 



dtp{x,t) 



^d. 



^{^V'{x) + !)-\d..g'-{x) 

1 7 



p{x,t) 
(5) 



-dx3{x,t) 



in Ito's sense where j{x,t) is the current and p{x,t) the 
probability distribution. 

If we first take the overdamped limit and afterwards 
assume Gaussian white noise, we end up with a Fokker- 
Planck equation in Stratonovich's sense. Both equations 
differ in a drift term which can be absorbed in an effective 
potential. Thus, the optimal current does not depend 



on the interpretation. The optimal potential is merely 
shifted by the drift term. In the following we use the 
Fokker-Planck equation in Ito's sense. 

We next introduce dimensionless units. Wc express 
energies in units of fcsTo where 



To = - / T{x) dx. 
^ Jo 



(6) 



By introducing a scaled length x = x/L and a scaled 
time t = t/tf) with ig = L'^^/ksTo, we rewrite the Fokker- 
Planck equation and identify the dimensionless quantities 



f{x) 



_T{x) 



To 



V{d 



_ V{a 



ksTo 



I 



fL 
ksTo' 



Kx, t) = — -— , p{x, t) = Lp{x, t). 

Kb To 



(7) 



For ease of notation we drop the hats in the following. 
The dimensionless potential V{x) and the dimensionless 
temperature profile T{x) are periodic, V{x -f- 1) = V{x)^ 
T{x -I- 1) = T{x). In the steady state, the current j is a 
constant and the Fokker-Planck equation reduces to 

J = [-V'{x) +f~ T'ixMx) - T{x)dMx)- (8) 

For a periodic temperature we solve this equation under 
the condition of a periodic p{x). Without loss of gener- 
ality, we chose T^(0) = which results in 



p{x) : 
with 



je 



1 



Tix) \l-e-<PWJo Jo 

(9) 



^^ ' ' T{x') 



(10) 



Using the normalization J„ p{x)dx = 1, we obtain the 
inverse current 



1 - e-^(i) Jo 



("') dx' 



1 e<^(^) 
Tix) 



dx 



e ^(^') dx' ' ' 



dx (11) 



/o \Jo / T{x) 

which has been derived previously in Refs. |20l.[33|. 

III. OPTIMIZING THE CURRENT 

The inverse current [Eq. (fTTj) ] depends on the shape 
of the potential V{x). Instead of optimizing the current 
directly with respect to the potential wc introduce 



B{x) = / e-^^^') dx' 

'0 



(12) 



and rewrite the inverse current in a more elegant way as 
a functional of B{x) 



r'[B,B']= f 

Jo 



with 



^ a~B{x) 
T{x)B'{x) 

^ B{1) 

i^B'iiy 



dx 



(13) 



(14) 



For the minimization of the functional ([T5)) . the function 
space is constrained by the periodicity of the potential 
which imposes a constraint on B{x). With the boundary 
condition V{1) = V{0) = we obtain from Eq. ^ the 
non-local constraint 



/ <l>'Tdx^f 
Jo 



(15) 



which by using Eq. (|12p can be transformed into the 
isoperimetric constraint 

»i 



T'lnB' dx = f + TlnB' 



(16) 



In order to minimize the inverse current [Eq. (J13p ] un- 
der the constraint [Eq. P^ ] we introduce the effective 
Lagrangian 



£[B{x),B'{x),x] 



a — B{x) 
T{x)B'{x) 



\T'{x)\nB'{x) (17) 



with a Lagrange multiplier A. 

For a unique solution to the corresponding Euler- 
Langrange equation we have to impose two boundary 
conditions. The first one, B{0) = 0, arises naturally from 
the definition of B{x) in Eq. (fT^ . One is tempted to use 
the condition -B'(O) = 1 as the second one, but this is not 
appropriate. In principle, we have to allow the derivative 
to jump at the boundaries because such jumps do not 
contribute to the integral of the Lagrangian. In this way 
the boundary condition fixes the value B'{0) = 1, but the 
value of linie^o -S'(0+|e|), which is the relevant boundary 
condition for the solution of the Euler-Langrange equa- 
tion, in fact, remains a free parameter. This feature has 
been discussed in detail for similar optimization problems 

113, mi. 

For the optimization we proceed in two steps. First 
we minimize the integral of the Lagrangian L C dx and 
then we adjust the remaining parameters to obtain the 
maximum current. The corresponding Euler-Langrange 
equation is given by 

2B'^T+{a-B){T'B'+2TB")+XT^B'{T"B'-T'B") = 

(18) 
with the boundary condition B{0) ~ 0. Changing vari- 
ables B{x) = — aexp[/(x)] + a leads to a second order 
differential equation for I(x) which is integrable. The 
solution 



/±(x) 



2dx' 



-XTT' ± ^\T{c -K AT) -h A2T2T'2 



(19) 



still depends on the Lagrange multiplier A and the new 
constant c. The solution of the Euler-Lagrange equation 



B{x) = — aexp[/±(x)] 
leads to the optimal inverse current 

dx 



(20) 



r'(A,c)=- 



1 



nx)iL{x) 

1 r 







4(|, + A)+A2r'2 



dx (21) 



which depends on the two parameters A and c but is 
independent of a. Note that in order to obtain a positive 
current, we chose the minus sign of the square root in 
Eq. dUl). 

In the next step, we optimize the inverse current [Eq. 
pTjl ] by adjusting the free parameters A and c. These 
parameters are not independent but related by the con- 
straint [Eq. HH)] 



n(A,c) = Tln|r 



+/- 



T'lnir I dx- 



I'_T dx = 0. 



In the appendix, we show that in general the optimization 
problem has a solution for c = 0. Hence, for a given 
temperature profile the remaining parameter A can be 
determined by the constraint [Eq. ((22)) ]. 

From Eqs. (fTUlll2[[^(I)) we derive the optimal potential 



V{x) ^T\n\r_ 



fx- 



I T'\n\I'_ 
Jo 



dx' 



I'_T dx' 

(23) 



which becomes the basis for the following case studies. 



IV. CASE STUDY I: SINUSOIDAL 
TEMPERATURE PROFILE 

In this section we will discuss a sinusoidal temperature 
profile 



T{x) = Asa\{2iTx) + l 



(24) 



with the amplitude Q < A < 1. The external force / is 
a second parameter. In the following we will study how 
the optimal potential and the current depend on these 
two parameters. 

The height of the potential is the essential blocking 
mechanism in the ratchet. The thermal "kicks" from 
the environment move the particle over this barrier. We 
expect the largest slope of the optimal potential roughly 
at the hottest point because there the fluctuations are 
strong enough to push the particle against a large force. 
In the colder regions, the potential should decrease. The 
optimal potential as determined numerically using Eq. 



indeed fulfills these expectations, see Fig. 2(a) 
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FIG. 2: (a) Optimal potential V{x) and T(x) = i sin(27ra;) + 1 for load / = -0.05. (b) Optimal potentials V(x) for T{x) = 
ylsin(27ra;) + 1 with different amplitudes A and load f — 0. Inset: Optimal current j versus A. (c) Optimal potentials V{x) 
for T{x) = Asm{2TTx) + 1 with different amplitudes A and load / — —0.05. Inset: Optimal current j versus A. (d) Optimal 
potentials V{x) for T{x) — i sin(27ra;) + 1 and for different loads /. Inset: Optimal current j{f) and power P{f) as a function 



of the force. The power exhibits a maximum Popt — 0.52 at /o 



opt 



-1.27. 



Optimal potential for different amplitudes of 
the temperature profile 



The optimal potential for different amplitudes A with 
zero external force / = is shown in Fig. 2(b) For 
amplitudes A ^ 1, the temperature in the colder area 



X > 0.5 goes to zero and the optimal potential becomes 
strongly asymmetric. In this low temperature area the 
thermal fluctuations are so weak that even a gently de- 
clining potential is sufficient to push the particle in one 
direction. The optimal current is proportional to the ab- 
solute value of the amplitude, see Fig. 2(b) This general 



scaling behavior is not limited to a sinusoidal tempera- 
ture profile which can be understood as follows. With 
the periodicity r'(0) = T'{1), the first term of the con- 
straint [Eq. dm] vanishes. With c = and / = 0, Eq. 



dm and Eq. dUl) become 



n(A) = / T' In AT' + ^/iX + \^T'^ 
Jo 



k-^ 



A2r'2 dx = Q (25) 



and 



J 



r\\)^ 



v/4A + A2T'2 dx. 



(26) 



The derivative of the temperature scales like T'{x) ex A. 
Choosing A oc A~^ , the constraint [Eq. d^5)) ] is indepen- 
dent of A. The current [Eq. ([^5]) ] then is a linear function 
of A 

Now wc consider the system with a finite external force, 
/ = —0.05. For temperature amplitudes A ^- 1, effec- 
tively corresponding to a large temperature difference, 



the external force can be neglected and the optimal po- 
tential looks like in the case without an external force 
see Fig 



2(c) 



For small A, the ratchet effect induced by 
the temperature difference is not strong enough against 
the external force. In this regime, the optimal potential 
has to prevent the particle from being dragged in the di- 
rection of the external force. This is achieved by blocking 



the particle with a larger barrier, see Fig. 2(c) 



B. Optimal potential for different external forces 

For stronger external forces, the potential has to com- 
pensate this dragging mechanism with a larger barrier 
in order to obtain a current in the direction opposite to 
the force. Thus, we expect larger potentials and smaller 
currents for stronger external forces which is confirmed 
by our calculations, see Fig. |2(d)[ 

We next calculate the dimensionless power output of 
the heat engine which is given by 



P = -fj- 



(27) 



The power output as a function of the force / is shown 
in the inset in Fig. |2(d)[ It exhibits a maximum at 
intermediate forces where the heat engine thus works in 
a maximum power regime. 



V. CASE STUDY II: PIECEWISE CONSTANT 
TEMPERATURE 

We now consider a piecewise constant temperature 



T{x) = 1 + AT - 2ATe{x - 1/2) 



(28) 



with a hot and a cold area with temperatures T^ — l+AT 
and Tc = 1 - AT, respectively, where < AT < 1. We 
then face discontinuities in the Fokker-Planck equation. 
In this section, we first analyze a continuous approxima- 
tion to the piecewise constant temperature. The optimal 
potential then has an complex shape with peaks. We 
compare these results with a numerical solution for the 
optimal potential. 

The continuous temperature profile 



VTTdsin(27r2;) 
1 (x) ~ ^^ + 1 



lJsm'^{2'!Tx) +d 



(29) 



interpolates between the extreme values d — > which 
corresponds to the piecewise constant temperature [Eq. 
(|28)) ] with AT = 1/2 and d ^ oo where the pr ofile b c- 

with A = 1/2, see Fig. 



3(a) 



comes sinusoidal [Eq. 

The optimal potential for different values of the pa- 
rameter d is shown in Fig. |3(b)[ For d <^ I, two peaks 
emerge in the optimal potential at the positions of the 
temperature discontinuities. In between these two peaks, 
the potential decreases linearly. The current diverges for 
d ^ 0, see inset of Fig. |3(b)[ 



A. Numerical solution for a piecewise constant 
temperature 

We next investigate this complex shape of the potential 
and the divergent current in more detail by solving the 
problem numerically on a discrete lattice. The goal is to 
minimize the inverse current [Eq. (jll|) ] for the piecewise 
constant temperature [Eq. (pS)) ] . The periodic boundary 
condition [Eq. (|15p ] for the potential transforms to the 
condition 



(1) = <^(l/2) 



2AT 
AT- 1 



/ 



1 -AT 



(30) 



considering that 0(0) = by definition [Eq. pO)) ] . For a 
numerical solution we discretize (f>{x) — *■ (j){xi) = (pi with 
i ~ 0, . . . ,N and search for a global minimum oi j~^{(t)i) 
in this N + 1-dimensional space with a simplex algorithm 
proposed by Nelder and Mead [39| . The boundary values 
00 and (/)Ar are given according to Eqs. pO l lSO l) by 



0, 



4>N^4>N/2 I ^y-^ 



1 -AT 



(31) 
(32) 



and 01, ... , (pN-i are varied to yield a maximum current. 
Using Eq. (fTU]) . we then calculate the potential V{x) 
from the optimal 0(x). 

The numerical solution for the optimal potential de- 
pends on the discretization which determines how large 
the gradient of the temperature and the potential can 
be, see Fig. 3(c) For finer discretization, the optimal 
potential shows larger gradients. As a consequence, we 
find that the current inum{N) diverges with increasing 
discretization N -^ oo, see Fig. |3(d)[ This is consi stent 
with the developing divergence visible in Fig. 3(b) 



B. Origin of the divergent current 

Due to the lattice discretization, the temperature pro- 
file can be considered to be linear with gradients with 
absolute value l/(2e), see Fig. U) We calculate the op- 
timal inverse current [Eq. (|2ip ] for such a temperature 
profile 

j~i(A, e) = VA2 + 16Ae2 + \/a(1 - 4e). (33) 

The Lagrange multiplier A and the discretization param- 
eter e are related by the corresponding constraint [Eq. 

(Ell 
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16e2 4e - 1 



VA 



+/ - 0. 

(34) 

In the appendix, we show with a perturbation method 
that for e — > 0, the constraint [Eq. ([M]) ] can be fulfilled 
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FIG. 3: (a) Continuous approximation T{x) [Eq. (|29p ] to a piecewise constant temperature for different parameters d. In the 
limit d -^ oo, T{x) approaches the sinusoidal profile T{x) = i sin(27ra::) + 1. (b) Optimal potentials V{x) for T{x) from Fig. 
3(a) with external force / = —0.05. Inset: current j versus parameter d. (c) Optimal potentials obtained from the numerical 
minimization on a discrete lattice for a piecewise constant temperature profile T{x) [ Eq. (|28p ] with AT — 1/2 and with an 
external force / — —0.05 for different discretization A'^. (d) For an external force / — —0.05, optimal current jnum as obtained 
from the numerics with discretization 2e — 1/N compared to the perturbative calculation developed in the appendix. 
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FIG. 4: Temperature profile on the lattice with gradients de- 
pending on e. 



if A ^ 0. For an approximated dependence A(e), we use 



an iterative method and obtain a sequence of functions 



A„+i(e) 



1 



(ln|A„(e)|-21n|2e|-l + /)2 



(35) 



with Aq = 1. For the first members we calculate the 
current [Eq. ([33]) ] as a function of e. These results are 
compared to the numerical current, where the discretiza- 



tion 1/N corresponds to 2e, see Fig. 3(d) The divergent 



behavior of the current obtained from perturbation the- 
ory is in good agreement with the divergence from the 
numerics. The discrepancy is on one hand due to the 
numerical integration on the lattice in Eq. pip and on 
the other hand due to the approximations implicit in the 
perturbation method. 

More insight into the origin of the divergent current 
can be gained from analyzing the probability distribution 
[Eq. ^] which can be calculated from the Fokker-Planck 
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FIG. 5: (a) Probability distributions for a rising peak with height [/ in a piecewise constant temperature T(x) [Eq. (|28|) ] with 
AT = 1/2 and load / — —0.05. (b) Sawtooth potential with superimposed peak and its corresponding probability distribution 
in a piecewise constant temperature T{x) [Eq. (|28|l ] with AT = 1/2 and load / — —0.05 for different peak heights. 



equation for a model potential. Guided by the optimal 
potential obtained from the numerics, see Fig. 



3(c) 



consider a sawtooth potential with a superimposed finite 
peak at the temperature discontinuity at x = 0.5. Note 
that we use this potential only as a case study aiming at 
a deeper understanding of the divergence of the current. 
The full optimal potential has a second peak at a; = 
and both the height of the peaks and the linear slopes be- 
tween the peaks diverge. A finite peak potential is not a 




FIG. 6; Peak as a limiting process. Each side is in one tem- 
perature region. 

(5-function, but can be considered as a result of a limiting 
process from a triangular shape, see Fig. [51 It is impor- 
tant that the peak arises at the discontinuity of the tem- 
perature profile, namely that the rising side of the peak 
is in one temperature area and the other side is in the 
other temperature area. Otherwise the peak would not 
contribute to the integrals in the current. Since the peak 
has infinitesimally small width, the forward rate is not 
decreased by a higher peak potential, contrary to a ris- 
ing barrier with finite width. A finite peak in a constant 
temperature area would not contribute to integrals and 
thus would have no effect on mean first passage times. 
In contrast, for a dichotomous temperature profile and a 
finite peak at the interface, mean first passage times in- 



volving a crossing of the peak are affected by its height. 
The backward mean first passage time increases exponen- 
tially with the peak height. In contrast and somewhat 
counter-intuitively, the forward mean first passage time 
decreases with increasing peak height due to the strong 
suppression of peak recrossings from the cold side. 

A finite peak with height U superimposed on a flat 
potential causes a depletion of the particles on the high 
temperature side and an accumulation on the other side. 
With increasing peak height U, this behavior saturates, 
see Fig. 5(a) In a sawtooth potential, the particles ac- 
cumulate in front of the barriers, see Fig. |5(b)[ If we 
combine both effects by superimposing a finite peak on 
a sawtooth potential, the depletion can compensate the 
accumulation, see Fig. 



5(b) 



The compensation of the accumulation is the main rea- 
son for the divergent current. For a sawtooth potential 
with large amplitude, the particles would accumulate in 
front of the barriers. The peaks counteract this effect 
and allow the particles to overcome the barriers. In be- 
tween the peaks, the particles (on average) follow the 
steep potential with large mean velocity. For finer dis- 
cretization, the peaks get larger and the potential in be- 
tween steeper corresponding to a stronger force which 
increases the mean velocity and thus the current. 

This behavior is only valid under the assumption of an 
overdamped dynamics where the particle instantaneously 
adjusts its velocity according to the Langevin equation 
((4]). In the underdamped case, the particle needs a cer- 
tain distance to reach the velocity corresponding to the 
local force. The overdamped limit is a good approxi- 
mation as long as the relaxation time of the momentum 
Tij = 771/7 is small compared to the (mean) time Tc it 
takes for the particle to cross a region with basically con- 
stant force. Considering the potential above with dis- 
cretization e, both the peaks and the linear slopes be- 



twccn the peaks should not become too large for the over- 
damped limit to be appropriate. The peaks constitute 
the largest slope in the optimal potential and therefore 
are crucial for the appropriateness of the overdamped de- 
scription. Considering a peak with height U and width 
2eL, the relaxation time tr should be small compared to 
Tc = eL/v" where v^ = \V'/^\ = U/{eLj) is the (mean) 
stationary velocity. Therefore, the overdamped limit is 
approriatc for 



m L^e^7 
— < —r^ or 

7 U 



m e 



(36) 



For a given value of r7i/(7^L^), the discretization e thus 
cannot become too small for the overdamped limit to be 
still valid. Since the divergence of the current occurs with 
decreasing values of e, the current presumably does not 
diverge in any realistic system, where the underlying dy- 
namics is underdamped. This question, however, is hard 
to decide conclusively, since the current for underdamped 
dynamics can only be determined numerically for a given 
potential. The optimization of the potential (with an in- 
finite number of degrees of freedom) is a computationally 
difficult task to be reserved to future work. For a colloidal 
particle of radius R ~ 1/xm in a temperature profile with 
To ~ SOOiC and periodicity L ~ 50/.im, we get 



7^L^ 



2- 10 



'11 



1 



ksTo 



(37) 



where we have used Stokes friction 7 = Giri-jR and the 
mass m = AirpR'^/S with viscosity t] ~ 10~^ Pa s and den- 
sity p ~ 10'^ kg/m'^. Even for the finest discretization e ~ 
5 • 10"'' used in our study, we have e'^/U ~ 3 • 10" VfcsTo 
and condition ([36]) is fulfilled. The overdamped descrip- 
tion thus is valid for the optimal potential for any real- 
istic (finite) temperature gradient. The (large) currents 
shown in Fig. 3(d) thus are in principle observable in ex- 



periments, provided the necessary temperature jump can 
be generated on the scale of 2e • L ~ 10"'^ • 50/im ~ 50nm. 
A genuine divergence of the current, on the other hand, 
may be prohibited by the onset of inertia effects. 



VI. CONCLUSIONS 

Using variational calculus, we have developed a 
method to calculate the optimal potential which maxi- 
mizes the current in ratchet heat engines for a continu- 
ous temperature profile. In the load free case, we have 
shown that the maximum current depends linearly on the 
amplitude of the temperature profile. 

In the case of a piecewise constant temperature, the 
current diverges for the optimal potential consisting of 
steep linear parts and peaks at the boundaries which in- 
duce a long range effect in the probability distribution. 
However, the underlying assumption of an overdamped 
dynamic is limited by the slopes of the optimal potential 
presumably resulting in a bounded current for a parti- 
cle with finite mass. In addition, a piecewise constant 



temperature is an artificial model rather than physically 
feasible. For any future nano- or microfluidic realization, 
the temperature gradients will be finite and thus the cur- 
rent. 

In principle, external potentials for a Brownian particle 
can be realized by laser traps. However, it might be very 
difficult to model the optimal potential in every detail. In 
particular for the dichotomous temperature profile, it is 
clear that finite peaks must be approximated by a barrier 
with a finite width in any experiment. In order to esti- 
mate the observable velocity, we consider a colloidal par- 
ticle with radius R trapped in the optimal potential for 
a sinusoidal temperature profile. In recent experiments, 
temperature gradients AT/L ~ 10^ K/m have been gen- 
erated [25. In the load free case, the optimal dimension- 

2A, see Fig 



2(b) 



less current is roughly j 

the scaled temperature amplitude AT/Tq 

we get a rough estimate for the stationary velocity 



where A is 
WithEq. © 
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100 nm 



SnrjRL R/nm s 



(38) 



under the assumption of Stokes friction with viscosity 
77 ~ 10"'^ Pa s. Although such a transport effect is small 
at a micrometer scale, future realisations at a nanometer 
scale will yield observable velocities. 

Our study can be extended in several directions. In 
principle, it would be interesting to calculate the effi- 
ciency at maximum power for our model and compare 
it to previous results where only the load force was op- 
timized. However, it is difficult to define efficiencies for 
continuous temperature profiles. Moreover, Hondou and 
Sekimoto pointed out in Ref. [4^ that heat transfer 
cannot be treated appropriately within the overdamped 
Langevin equation. For a similar model of a ratchet heat 
engine, a heuristic argument was used to prop ose a po- 
tential which leads to a large Peclet number [4l| . It would 
be interesting to see whether our approach can be used to 
calculate the optimal potential with respect to a maximal 
Peclet number. 

While ratchet heat engines have not been realized ex- 
perimentally yet, the recent successful generation of large 
temperature gradients [23, [2J, [2^ may facilitate the con- 
struction of microscopic heat engines. By using our re- 
sults, such nanomachines may then be tuned to produce 
maximum power. 
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VII. APPENDIX 



constraint [Eq. ((22 



A. Optimization of the current with respect to the 
remaining free parameters A and c 

With I-{x) from Eq. (|19p . the inverse current [Eq. 
(EDI reads 



^ Jo 
with the constraint [Eq. 

XT' 



] for A and c 



dx (39) 



n(A, c) = / T' In 
Jq 

r-l 



4(- + A)+A2T'2 



dx 



-2dx 



AT' + ^4(f + A) + A^T'^ 



+ / = 0. (40) 



Note that the first term in Eq. ([22]) vanishes for r'(0) = 
T'(l) and r(0) = r(l). By introducing 



fc(A,c, ^)=j ^(A, c) -pin(A, c) 



(41) 



with Lagrange multipher fi we formulate a minimization 
problem under a constraint. For the optimal parameters 
c, A, fi which minimize the optimal inverse current [Eq. 
([3911 ] and fulfiU the constraint [Eq. (gO])], |f , || and 
n(A, c) have to vanish. These equations cannot be solved 
analytically. Nevertheless for c = we have 



dk 
'd'c 



(A-/.)/' 
Jo 



dx 



c=o Jo Ar(a;)V4A + A2T'2(x) 

which vanishes for X = fi. The partial derivative 



dk 
dX 



A— /j,,c— 



T'(a;) dx 



(42) 



(43) 



also vanishes for a periodic temperature profile. Thereby, 
we reduce the problem to a one dimensional root search 
of n(A, c)|c=o which can easily be done numerically. 

Thus, we find one solution but we cannot exclude rig- 
orously that there exist other solutions for this minimiza- 
tion problem. In our case studies, we find (by a compar- 
ison to numerical optimization) that the solution with 
c = indeed is a global minimum. This strongly sug- 
gests that it generally is the relevant solution. 



B. Perturbation theory for the divergent current 

We approximate the dependence A(e) for a tempera- 
ture profile as shown in Fig. 0] in the limit e — > 0. The 



In 



. A 1 /A2 ^^ 
8e2 4e V 4e2 



-\ I 



'''%'-^H - 



VA 



relates A with e. We define h{e) by 
A(e) EE Ae^h{e) 
and from Eq. (|^^ follows 
h 1 



Vh\n\l+- + -Vh^ + Ahl-Vh + 4+2+ fVh : 



2e' 



(44) 
(45) 

(46) 



For e — > 0, the right hand side of Eq. (|46p diverges and 
thus the left hand side must also diverge, yielding ft. — > oo. 
We consider the leading terms in Eq. (|15)) 



^/ft(ln|ft|-l + /) = l 



and rewrite it in a self consistent equation 

1 



4e2(ln|ft|-l + /)2- 
We use Eq. (|^5)) to obtain the corresponding 

1 



X{e) 



(ln|^|-l + /)2- 



(47) 



(48) 



(49) 



In perturbation theory, iteration is a standard method 
for algebraic equations [42| . In this way we search for 
a fixed point which is a solution for the equation. We 
iterate Eq. ((49|) 



A„+i(e) 



1 



(ln|A„(e)|-21n|2e|-l + /)2 
and choose 

Ao = l 
for a good convergence. More iterations lead to 

Ai(e)= ^ 



(21n|2e| + l-/)2 



A2(e) 



(2 In |2 In |2e| + 1 - /[ + 2 In |2e| + 1 - /)2 



(50) 



(51) 



(52) 



(53) 



and so on. Note that this sequence converges slowly but 
still gives an idea how A(e) behaves. For each A we ob- 
tain a current from Eq. ()33p . The convergence of this 
sequence of currents is shown in Fig. 3(d) 
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